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Abstract 


The HZETRN deterministic radiation transport code is one of several tools de- 
veloped to analyze the effects of harmful galactic cosmic rays (GCR) and solar 
particle events (SPE) on mission planning, astronaut shielding and instru- 
mentation. This paper is a comparison study involving the two Monte Carlo 
transport codes, HETC - HEDS and FLUKA, and the deterministic transport 
code, HZETRN. Each code is used to transport ions from the 1977 solar min- 
imum GCR spectrum impinging upon a 20 g/cm 2 Aluminum slab followed by 
a 30 gj cm 2 water slab. This research is part of a systematic effort of verifica- 
tion and validation to quantify the accuracy of HZETRN and determine areas 
where it can be improved. Comparisons of dose and dose equivalent values at 
various depths in the water slab are presented in this report. This is followed 
by a comparison of the proton fluxes, and the forward, backward and total neu- 
tron fluxes at various depths in the water slab. Comparisons of the secondary 
light ion 2 H, Z H , 3 He and 4 He fluxes are also examined. 


1 Introduction 

As part of the NASA strategic plan for the human exploration of space and the desire 
to provide radiation protection for both astronauts and instrumentation, there is ongoing 
development, testing, verification and validation of the deterministic radiation transport 
code HZETRN (High charge (Z) and Energy TRaNsport), which was developed at Langley 
Research Center over the last twenty - five years by Wilson and co - workers [1]. The 
present paper, which focuses on galactic cosmic rays (GCR), is an extension of reference 
[2], which involved a verification of HZETRN for solar particle events (SPE). 

HZETRN is a deterministic computer code which has been used for radiation analysis 
under a variety of shielding conditions in analyzing solar particle events, galactic cosmic 
rays, and low earth orbit (LEO) trapped proton environments. While the code has en- 
dured several rounds of verification and validation in these environments [1, 3, 4], most 
of the comparisons have focused on integrated quantities, such as dose (D) or dose equiv- 
alent (H), and individual ion fluences were not examined in detail. Though dose or dose 
equivalent were generally viewed as sufficient tests for evaluating code accuracy, recent 
interest in fluence based approaches to radiation risk assessment demand a higher degree 
of accuracy [5]. This, along with recent improvements to some of the underlying trans- 
port models and numerical procedures within HZETRN [6, 7], provides an opportunity 
for another round of verification and validation benchmarks. 

In this paper, the 1977 solar minimum GCR spectrum is used to compare HZETRN to 
the Monte Carlo codes HETC - HEDS (High Energy Transport Code - Human Exploration 
and Development of Space) [8, 9] and FLUKA (FLUctuating KAskade) [10, 11]. Models 
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of the free space GCR environment developed in the past two decades have provided the 
most realistic description of the interaction of incoming GCR from outside the heliosphere 
with solar activity. For the 1977 GCR minimum, the model of reference [12], was used 
as GCR input into the HZETRN, HETC and FLUKA code systems. This GCR model 
is based on fitting the existing balloon and satellite measured differential energy spectra 
from 1954 - 1992 to the stationary Fokker - Planck equation to estimate the appropriate 
diffusion coefficient. The implementation of this model accurately accounts for the solar 
modulation of hydrogen through nickel (H - Ni) by propagating the local interplanetary 
spectrum of each element through the heliosphere. This is accomplished by numerically 
solving the Fokker - Planck diffusion, convection and energy loss boundary value problem. 
The model provides a single value of the deceleration parameter, 0(£), describing the level 
of solar cycle modulation and determines the GCR differential energy spectra for all of 
the elements at a given radial distance from the sun. 

The benchmark problems solved by the three computer codes consist of a 20 g/cm 2 slab 
of Aluminum shield followed by a 30 g/cm 2 water slab target when exposed to selected 
ions from the 1977 solar minimum GCR spectrum. The thickness of the Aluminum 
slab is representative of some spacecraft thicknesses, and the water slab thickness of 
30 g/cm 2 was selected because thicknesses between 0 and 30 g/cm 2 can be used, with 
appropriate interpolation methods, to determine body organ exposures. Due to their 
prominence in the GCR spectrum, the fluxes associated with Hydrogen, Helium, Oxygen, 
Carbon, Magnesium, Silicon and Iron were selected as individual boundary conditions to 
the Aluminum slab. 

The computational geometry is illustrated in Figure 1. Figure 2 illustrates the magni- 
tude of the flux boundary terms from the 1977 solar minimum GCR spectrum associated 
with the various ions of interest. In this study, each code is used to calculate specific 
radiation quantities at various depths in the water slab. In particular, each code provided 
the dose rate (cGy/day) and dose equivalent rate (cSv/day) at various depths in the wa- 
ter slab. Here, dose is the mean energy imparted (absorbed) per unit mass, and dose 
equivalent is a measure of biological effectiveness. This is followed by a comparison of the 
proton and forward, backward and total neutron fluxes. Comparison graphs for light ion 
flux calculations associated with 2 H, 3 H, 3 He and 4 He ions are examined. Also included 
are comparison graphs of selected fluxes associated with heavy ion production. 

The benefits of improving HZETRN, when two approaches (HETC, FLUKA) already 
exist, is that HZETRN is able to run much faster than these other Monte Carlo codes. 
Therefore, HZETRN is superior in terms of engineering and other trade studies. A de- 
scription of HZETRN, plus the recent improvements involving neutron and light ion trans- 
port, has been fully discussed in reference [2] and will not be repeated here. These recent 
improvements to HZETRN are used extensively in the present work. 
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2 Monte Carlo codes 

The term Monte Carlo, is used to describe a random walk algorithm that simulates an 
event and then performs a statistical analysis of the results. Monte Carlo techniques are 
typically constructed whenever there is a high dimensional integral equation that cannot 
be solved either analytically or numerically with deterministic methods. To reduce the 
computing time involved in Monte Carlo techniques, variance reduction is often employed. 
A general description of Monte Carlo methods can be found in reference [13]. 

In this paper, results from the two Monte Carlo codes, HETC - HEDS [8, 9] and 
FLUKA [10, 11], are used to compare with the results calculated by the deterministic 
code HZETRN. The following is a brief description of HETC - HEDS and FLUKA codes 
used in this comparison study. 

2.1 HETC - HEDS 

The HETC - HEDS computer code is a Monte Carlo based solution method designed 
specifically for solving space radiation problems [8, 9] associated with secondary particle 
fields produced by space radiation interacting with shielding and equipment. It is a 
three dimensional generalized radiation transport code capable of handling and analyzing 
radiation fields which affect critical body organs of astronauts such as bone marrow and 
the central nervous system. HETC - HEDS can be applied to a wide range of particle 
species and energies. It contains a heavy ion collision event generator capable of tracking 
nuclear interactions and performing statistical analysis of the data. It simulates particle 
interactions by using a pseudo random number generator, along with the appropriate 
physics, to follow trajectories of primary particles and all secondary particles involved 
in the nuclear collision of galactic cosmic rays and solar event particles interacting with 
shielding material, biological organisms and electronic equipment. The geometry input is 
simple combinatorial volumes. The cross sections are generated internally, as part of the 
software. The HETC - HEDS computer code employs all particles of interest for space 
radiation. In particular, HETC - HEDS considers interactions of protons, neutrons, 7r + , 
7T _ , /i + , light ions and heavy ions. These particles can be arbitrarily assigned position, 
angle and energy throughout a spatial boundary. This Monte Carlo code follows each 
particle in a cascade until it undergoes a nuclear collision, absorption, decay or escape 
from the spatial boundary, or is eliminated by crossing a domain variable cutoff. The 
nuclear reactions and processes are accounted for by using appropriate physical models 
to handle such things as energy loss, range straggling, Coulomb scattering, etc. Both 
elastic and nonelastic collisions are considered using energy and nucleon conservation 
principles. A user written part of the code generates a source which allows for a simple 
form of biasing. This Monte Carlo package has no capability to determine when the 
solution method has converged on the answer sought within some confidence interval. It 
does provide the means for predicting the interaction product yields, production angles 
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and energies using nuclear models for transport processes. It has been extensively used 
for code verification, benchmarking and testing against available laboratory beam data 
associated with energetic heavy ions. 

2.2 FLUKA 

The FLUKA computer code is a general purpose Monte Carlo computer program used 
for calculating particle transport and interaction with various materials. It has the ability 
to transport and interact all elementary hadrons, light and heavy ions, and electrons and 
photons over an energy range which extends up to 10 4 TeV for all particles, and down to 
thermal energies for neutrons [11]. The code has built-in capabilities for scoring particle 
fluences, yields, and energy deposition over arbitrary three dimensional meshes, both on 
an event - by - event basis or averaged over a large number of histories. This code has 
been extensively benchmarked against available accelerator and cosmic ray experimental 
data, at beam energies as low as a few MeV and as large as cosmic ray energies. 

The code has special add - on modules for the generation of GCR spectra. The 
spectra are based on a blend of a model originally developed at NASA and extensions 
and modifications aimed at achieving the best fit of the most recent available experimental 
data, in particular those of the AMS (Alpha Magnetic Spectrometer) and BESS (Balloon 
- borne Experiment with a Superconducting Spectrometer) experiments. These spectra 
can be modulated within FLUKA according to an arbitrary solar activity modulation 
parameter, or for dates in the past, using the actual solar activity as measured by ground 
based neutron counters. Three different Earth magnetic field descriptions of increasing 
complexity can be used. Geomagnetic cut - offs can be input or calculated. Spectra 
representative of some of the largest solar particle events are also built into the code. 

In addition to galactic cosmic rays and solar particle events, this code has a wide range 
of particle physics application areas where it can be employed. It has been demonstrated 
that the code can accurately simulate nuclear interactions and events with great accuracy 
in several fields, including particle physics, dosimetry, accelerator driven nuclear systems, 
detector characterization, aircraft crew dosimetry, shielding, and hadron therapy. The 
code can simulate interactions involving over sixty different particle types and can han- 
dle complex geometries. The geometry input has two levels: the basic one, consisting 
of combinatorial volumes, and a second level, the so called “lattice” level, which allows 
repetition and arbitrary spatial placement of the complex object defined at the first level. 
It also allows for a voxel (a volume element, representing a value on a regular grid in 
three dimensional space) input stream which can be combined with a standard combina- 
torial geometry input, allowing, for example, to embed a detailed voxel description of a 
human being, such as derived from a computed tomography (CT) scan, within an arbi- 
trary spacecraft shape. The code can account for magnetic fields of arbitrary complexity. 
The nuclear models and associated cross sections are hard - coded into the software for 
all particles and energies, with the exception of neutrons below 20 MeV, where a 260 
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neutron group library obtained from standard international evaluated data files is used. 
The FLUKA computer code contains many variance reduction algorithms and covers all 
particles of interest for space radiation research. For most applications, the FLUKA pack- 
age requires no additional programming. The code provides a large variety of statistical 
techniques for analyzing nuclear interactions, and it can handle only one source particle 
type at a time without writing a special subroutine to perform many particle GCR, SPE, 
and LEO sources. The code can be used in either a biased or unbiased mode. The offline 
analysis tools provided with the code enable one to perform simple statistical analysis of 
the results. 

The code comes with a powerful Graphical User Interface, FLAIR, which allows an 
easy setup of the input stream and the geometry. It also provides user friendly tools for 
running the code and analyzing and plotting the results. 


3 Analysis of dose and dose equivalent curves 

Dose and dose equivalent values for various depths in the water slab are presented in 
Figures 3-9. These curves show some large differences in the dose and dose equivalent 
calculations generated by the three transport codes. For a quantitative assessment of the 
dose and dose equivalent curves, we select a root mean square (RMS) analysis of the data. 
The root mean square between data sets {!*} and {yi} for i — 1, . . . N is given by 


RMS = 


- y ,) 2 

i= 1 


(i) 


where N is the number of ordinates in the data sets with the abscissa values being the 
same for each set. A RMS comparison between the HZETRN, FLUKA, and HETC - 
HEDS data sets for dose are presented in Table 1. A similar root mean square difference 
calculation is performed on the dose equivalent data set, resulting in the values given in 
Table 2. 

One reason for the discrepancies in the values obtained for the dose and dose equiv- 
alent curves is due to the different ways that they are calculated by the three codes. 
Another reason for differences is due to differences in the physics models used for cross 
section calculations. Another difference is due to different quality factors that are being 
used in the calculations for dose equivalent. The major differences in the dose and dose 
equivalent calculations at various depths in the water slab by the three codes can be better 
understood by an examination of the methods used for calculation of the dose and dose 
equivalent values. Note that the following description of dose and dose equivalent calcu- 
lations, employed by the three computer codes also applies to the previous comparison 
study [2]. 
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Table 1: Dose rate (Gy/day) RMS differences between codes. 


Element 

HZETRN, HETC 

HZETRN, FLUKA 

HETC, FLUKA 

H (Figure 3) 

0.0115 

0.0119 

0.0026 

He (Figure 4) 

0.0003 

0.0088 

0.0089 

C (Figure 5) 

0.0003 

0.0002 

0.0002 

0 (Figure 6) 

0.0003 

0.0002 

0.0003 

Mg (Figure 7) 

0.0001 

0.0022 

0.0020 

Si (Figure 8) 

0.0001 

0.0002 

0.0001 

Fe (Figure 9) 

0.0001 

0.0004 

0.0003 


Table 2: Dose Equivalent rate (Sv/day) RMS differences between codes. 


Element 

HZETRN, HETC 

HZETRN, FLUKA 

HETC, FLUKA 

H (Figure 3) 

0.0112 

0.0124 

0.0212 

He (Figure 4) 

0.0049 

0.0155 

0.0112 

C (Figure 5) 

0.0010 

0.0009 

0.0019 

O (Figure 6) 

0.0011 

0.0020 

0.0030 

Mg (Figure 7) 

0.0003 

0.0121 

0.0119 

Si (Figure 8) 

0.0001 

0.0009 

0.0010 

Fe (Figure 9) 

0.0005 

0.0037 

0.0043 


3.1 FLUKA calculation of dose and dose equivalent 

FLUKA calculates the dose by scoring energy deposition events by all transported 
particles in a given volume and dividing by the mass. The dose equivalent is not a physical 
quantity, and consequently its calculation can be done in several ways. The FLUKA code 
does not compute dose equivalent directly. The three ways that dose equivalent related 
quantities are estimated by FLUKA are now described. The FLUKA dose and dose 
equivalent calculations in this paper were calculated using method 3 below. 

3.1.1 Method 1 

FLUKA calculates ambient dose equivalent by folding particle fluences (in volumes or 
on a surface) with appropriate fluence - to - dose conversion coefficients [14, 15]. When 
fluences on a surface are used, they must be properly computed by counting the particles 
crossing the surface weighted by unity divided by the cosine of the incident angle. 
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3.1.2 Method 2 

This method involves the use of KERMA factors as defined by NCRP [16]. KERMA 
is a acronym for Kinetic Energy Released per unit MAss or Kinetic Energy Released in 
MAterial. Therefore KERMA is similar to Dose and has the same units of Gray (Gy). A 
KERMA factor represents the ratio of the summation of the initial kinetic energies of all 
the charged particles which are liberated by uncharged particles, per unit mass of specified 
material. KERMA can be calculated for any specified material, at a point in free space, or 
at a point in an absorbing medium. Calculation of dose equivalent involves the weighting 
of each energy deposition event by a proper quality factor, defined by Q(L) : given by 
the ICRP60 curve [17]. FLUKA has technical issues implementing this weighting with 
neutrons below 20 MeV where average KERMA, rather than recoiled spectra, is usually 
employed. A special treatment for this type of calculation has been devised and is suitable, 
but it is limited to dose equivalent calculations as defined by ICRP 60 [17] and must be 
implemented in a user routine. 

3.1.3 Method 3 

This method uses the weighting of charged particle fluences, in a volume or on a 
surface, by the quality factor corresponding to their linear energy transfer (LET). This 
approach again gives the ICRP 60 kind of dose equivalent [17], because the code has 
estimators for the LET spectra. However, this method implicitly neglects the neutron 
and partial photon contributions, since neutron recoils below 20 MeV are not explicitly 
generated (apart from recoil protons on H). Note also that recoils are often so low in energy 
as to be either below transport cut off or, given the short range, are hard to sample in 
a statistically meaningful way. Again, it is mandatory for a correct fluence estimation at 
surfaces to use inverse cosine weighted scoring. The FLUKA dose and dose equivalent 
calculations in this paper were calculated using this method. 

3.2 HETC - HEDS calculation of dose and dose equivalent 

The computer code HETC - HEDS employs other codes such as MORSE [18], EGS 
[19] and MCNP [20]. The histories generated by HETC - HEDS are stored in output 
files. Transport of neutrons below 20 MeV and electrons is not carried out by HETC 
- HEDS. Instead, the output file of neutrons having energies of 20 MeV and below is 
fed as input into MORSE or MCNP to be transported by those codes. Electrons and 
photons from the HETC - HEDS output file are fed into EGS for transport. All primary 
and secondary charged particles, except electrons, are transported by HETC - HEDS and 
have their histories contained in the HETC - HEDS output file. The energy depositions 
from HETC - HEDS and EGS are straightforward. They are calculated by examining the 
energy of the particle at one position and then at another, usually defined by boundary 
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crossings, collisions, or end of range. The energy deposited is then the difference in the 
energies at the two sites. This is a standard Monte Carlo approach. Since the particle 
type and energy are known, the LET quality factor Q, [17], can be applied to convert 
from Gray to Sievert. For the low energy (< 20 MeV) neutrons, the flux is calculated as 
follows. The flux - to - energy deposited KERMA factors are integrated with the flux to 
obtain Gray. The neutron energy dependent quality factors are used with the KERMA 
factors to convert the Gray to Sievert. The KERMA factors used by HETC are those for 
water. For the proton source data, where the neutron dose is important, the outdated 
neutron energy dependent quality factors were used instead of the most recent ones, which 
is another source of error for the differences in dose equivalent. For the heavy ion data, 
the low energy neutron doses are of less importance. It should be noted that new Sievert 
data, using the more recent quality factors, was provided for the solar flare results [2], 
and the revised data compared favorably with other calculations. 


3.3 HZETRN calculation of dose and dose equivalent 

In HZETRN, absorbed dose D, due to energy deposition at a given location x by all 
particles, is calculated according to 

COO 

J>(*) = £/ S j (E)^ j (x,E)dE (2) 

j J0 

where (f)j(x, E ) is the fluence or flux of type j ion at depth x with kinetic energy A, and 
Sj(E) is the stopping power of type j ion with energy E. 

For human tissue exposure estimation, HZETRN calculates the dose equivalent H in 
terms of the product of dose D and quality factor Q, where the ICRP 60 quality factors 
[17] are based upon cancer risks. These factors relate the biological damage incurred from 
any ionizing radiation to the damage produced by the reference soft X - ray exposure. 
By definition, Q is a function of LET of the radiation field, which depends on both the 
ion type and its energy. The exact functional form of Q is prescribed in the process of 
setting radiation guidelines as defined by the ICRP publication 60 [17], and expressed 
below; hence, Q is a legislated quantity rather than the result of a particular measure- 
ment. However, the dependence of Q(L ) on the value L = LET, is intended to reflect 
a judgment related to the dependence of relative biological effectiveness (RBE) on LET 
at the anticipated level of exposure. In calculation of the dose equivalent, the HZETRN 
code uses the linear energy transfer quality factor Q(L ), as defined in reference [17] 


Q{L) = 


1 , 

0.32L - 2.2, 
300 /VL, 


L < 10keV//wn 
10 < L < 100 keV / //m 
L > lOOkeV/Twn 


( 3 ) 
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The graph of Q(L ), is illustrated in Figure 10. Note the difference in the quality factors 
for the ICRP 60 [17] and ICRP 26 [21] in the region 10 < L < 1000 kev///m, where the 
ICRP 26 values are given by [21] 


f 

Q{L) = \ 13.71og(l + L/466), 

20 , 


L < 1 keV//im 
1 < L < 20 keV / (Km 
L > 20 keV / (Km 


( 4 ) 


The value of dose equivalent is computed at a given location x from all particles present 
according to the relation 


roc 

H(x) = £ / Q s {E)Si{E)*fa E ) dE ( 5 ) 

j J0 

4 Comments on dose and dose equivalent 

The HZETRN, FLUKA and HETC - HEDS transport codes use different ICRP stan- 
dards for calculating dose and dose equivalent as illustrated in Figure 10 [13, 16, 22], 
HZETRN does not transport certain particles such as pions, muons, positrons, electrons 
or photons. These particles are used in calculating dose and dose equivalent by both 
HETC - HEDS and FLUKA. The contribution of these particles to dose and dose equiva- 
lent values can be significant. For low energy neutrons, the HETC - HEDS code uses the 
NCRP 132 and 142 [16, 22], conversion factors in combination with the DABL69 KERMA 
factors [23]. The HZETRN code does not use KERMA factors in the calculation of dose 
and dose equivalent. HZETRN transports light target fragments explicitly, so all of the 
neutron dose and dose equivalent is coming from target fragmentations and recoils with 
charge greater than 2. HZETRN uses a parameterization to do this. (No documentation 
exists for this parameterization. Future work must correct this.) Also, note that the 
FLUKA code neglects the contribution from low energy (< 20 MeV) neutrons altogether. 
The number and type of particles used and the algorithms employed by the three codes for 
calculating the dose and dose equivalent values are the probable causes for the differences 
in the dose and dose equivalent curves [13, 16, 22, 23]. 

5 Flux calculations 

A quantitative measure of “goodness of fit” for each of the flux comparison curves 
found in Figures 11 - 65 is constructed for comparison purposes. At any fixed energy 
value where the flux values from HZETRN, HETC - HEDS, and FLUKA are well defined, 
one can calculate the absolute value of the percent errors between HZETRN and HETC 
- HEDS, as well as the absolute percent error between HZETRN and FLUKA. These two 
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percent error absolute values are then added to get a sum of the percent errors at the 
selected energy value. We define this as the percent error norm at the selected energy 
value £ , and use the notation || E ||% to denote this sum of errors at energy E. Each of 
Figures 11 - 61 is examined and for each set of curves, three energy values are selected. 
The first energy value £1 is selected in the range 1 < £i < 10 MeV, the second energy 
value £2 is selected in the range 10 < £2 < 100 MeV, and the third energy value £3 is 
selected in the range £3 > 100 MeV. The values {£1, £2, £3} are selected to maximize 


£ = Int 




% + || £2 ||% + || -E3 


( 6 ) 


where the integer calculated for £ can be interpreted as a weighted average worst case 
percent error value assigned to a set of graphs. These F value error calculations are sum- 
marized in Tables 3-4. The blank portions in Tables 3-4 indicate that flux calculations 
were not provided for comparisons by either HETC - HEDS or FLUKA. 

Compare the £ values in Tables 3-4 with an examination of the set of curves in 
the corresponding figure number and note that for £ values less than 100 the agreement 
between the curves is quite good. For £ values satisfying 100 < £ < 1000, there are 
regions within the graphs where the three curves differ and for £ values greater than 
1000, there are major differences in the three sets of curves. The major differences are 
possibly the result of using inaccurate cross sections in performing the calculations. This 
hypothesis can be investigated and tested at a later date. 

Table 5 gives the average £ values for the sets of flux graphs corresponding to the 
forward neutrons, backward neutrons, total neutrons, protons, 2 H, 3 H, 3 He and 4 He. 


6 Additional graphs 

The light ion contribution is one of the major contributors to the dose and dose equiv- 
alent calculations. It is also known that there are some charge removal processes resulting 
from the collision of heavy ions which can also contribute to the dose and dose equivalent 
calculations. Figures 66 - 69 represent nuclear fragmentation effects in transport processes 
using the spectra of GCR Carbon, Oxygen, Magnesium and Iron. For example, Figure 
66 illustrates transmitted Carbon, Boron and Beryllium flux versus energy at the Alu- 
minum - water interface from GCR Carbon ions incident upon the 20 g/cm 2 Aluminum 
slab. Similar graphs are given in the Figures 67 - 69. The £ - norm for these set of 
curves are given in Table 6. Clearly, Beryllium production from Carbon needs the most 
improvement. 
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7 Conclusions 


The HZETRN space radiation code has been developed to analyze the effects of galac- 
tic cosmic rays and solar particle events on mission planning and shielding for astronauts 
and instrumentation. We have provided a comparison study involving two Monte Carlo 
transport codes, HETC - HEDS and FLUKA and the deterministic code, HZETRN. Each 
code is used to transport an ion from the 1977 solar minimum GCR spectrum impinging 
upon a 20 g/cm 2 aluminum slab followed by a 30 g/cm 2 water slab. We have compared 
dose, dose equivalent and fluence of a variety of particles at various depths in water. 

A parameter F has been calculated that can be interpreted as a weighted average 
worst case percent error value, and this value is assigned to a set of graphs depicting the 
transport code results. The large F value associated with a set of flux graphs indicates 
that there is significant difference between the HZETRN results and the HETC - HEDS 
and FLUKA results. A large F value will always occur whenever there is a significant 
spread between the curves HETC - HEDS and FLUKA. There are many instances where 
the curves for HETC - HEDS and FLUKA do not agree. This is possibly due to the 
cross sections used within these codes. Whenever HETC - HEDS and FLUKA disagree, 
further investigation and verification of results needs to be performed. The significantly 
high F values in Table 5 indicate that HETC - HEDS and FLUKA have large differences 
with HZETRN. Cross section calculations are a cause of error that need to be examined 
by updating cross sections models in all three computer codes. Another source of error is 
possibly due to how the light ion flux calculations are performed within HZETRN. These 
calculation methods need to be investigated in more detail. 

Figures 64 - 65 are magnified views of selected data sets illustrating that there are 
still major differences in the flux calculations as performed by HZETRN, HETC - HEDS 
and FLUKA, which in turn lead to large F values. These differences are perhaps due to 
the cross section values used during calculations. The magnified views are presented over 
the energy regions 1 < E < 500 MeV, because it is in this region that the quality factors 
have their greatest rate of change. These differences in values are another possible reason 
why the dose and dose equivalent calculations are not in agreement. 

Future research into a comparison of the results from the HZETRN, HETC-HEDS 
and FLUKA transport codes for the same shield/target configuration, will aid in trying 
to isolate causes for the differences between the codes found herein. Results from this 
study and future comparison studies can provide insight into the HZETRN code operation 
and bring the code to a high standard of performance so that design engineers can have 
increased confidence in its abilities to analyze particle transport associated with a variety 
of radiation environments. 
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Table 3: F values and quantitative comparison of flux graphs. Abbreviations are as follows: forward neutron (f - n), 
backward neutron (b - n), total neutron (t - n), proton (p). The decimal number after the figure number refers to 
the value of the water depth in g/cm 2 . For example, 11.0 refers to Figure 11 at a water depth of 0 g/cm 2 , which is 
the upper left panel of Figure 11. 
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30.30 



46 

32.30 

32 

33.30 


34.30 
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35.30 

1560 

36.30 

105 

37.30 

52 

30.20 
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31.20 

53 
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15 

33.20 
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34.20 
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35.20 

1745 

36.20 
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67 

30.10 
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31.10 

99 

32.10 

29 

33.10 

406 

34.10 

2704 

35.10 

2243 

36.10 

126 

37.10 
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31.0 

84 
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41 
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75 
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64 

24.30 

H 

H 
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26.30 
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27.30 

1698 

28.30 

242 

29.30 

77 

22.20 

422 
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37 

24.20 

ZO 

25.20 

311 

26.20 

2938 

27.20 

2679 

28.20 

316 

29.20 

87 

22.10 

199 

23.10 

47 

24.10 

oo 

25.10 

291 

26.10 
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27.10 

3892 

28.10 

539 

29.10 

91 

22.0 

146 

23.0 

47 

24.0 

00 

25.0 

372 

26.0 
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27.0 

14382 

28.0 
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29.0 
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72 

16.30 



101 

18.30 

99 

19.30 

471 

20.30 
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99 
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78 
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51 

19.20 
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65 
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17.10 

76 

18.10 

48 

19.10 

1227 

20.10 
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21.10 

135 

16.0 

66 

17.0 

98 

18.0 

55 

19.0 

1555 

20.0 
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21.0 
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50 

11.30 



59 

13.30 

31 

14.30 

2085 

15.30 
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46 
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29 

14.20 
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49 
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34 
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Table 4: Table 3 continued. 


Iron 

29 

54.30 



43 

56.30 

15 

57.30 

421 

58.30 

646 

59.30 

1096 

60.30 

128 

61.30 

34 

54.20 

953 

55.20 

40 

56.20 

10 

57.20 

480 

58.20 

798 

59.20 

1586 

60.20 

144 

61.20 

39 

54.10 

825 

55.10 

42 

56.10 

16 

57.10 

556 

58.10 

1072 

59.10 

2019 

60.10 

144 

61.10 

48 

54.0 

1003 

55.0 

47 

56.0 

16 

57.0 

657 

58.0 

2843 

59.0 

4570 

0'09 

165 

019 

Silicon 

47 

46.30 



47 

48.30 

26 

49.30 

391 

50.30 
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51.30 
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52.30 
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53.30 

57 
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45 
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LO 
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427 
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72 
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50 
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22 

49.10 
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2076 
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57 
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23 

49.0 
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51.0 
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52.0 
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39 
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43 
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28 
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86 
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45 
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46 

40.20 

24 

41.20 

373 
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43.20 
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44.20 

97 
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54 

38.10 
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55 

40.10 

22 

41.10 
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42.10 
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44.10 
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45.10 

95 

38.0 
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77 
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29 
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Table 5: Average F values. 


Forward Neutrons 

64 

Backward Neutrons 

753 

Total Neutrons 

56 

Protons 

26 

2 H 

630 

3 H 

3193 

3 He 

3001 

4 He 

310 


Table 6: F values heavy ions 


Figure 66, Incident ion is Carbon 

Carbon 

Boron 

Beryllium 

59 

66 

5850 

Figure 67, Incident ion is Oxygen 


Oxygen 

74 

Nitrogen 

106 

Carbon 

53 

Figure 68, Incident ion is Magnesium 


Magnesium 

78 

Sodium 

162 

Neon 

103 

Figure 69, Incident ion is Iron 


Iron 

215 

Magnesium 

241 

Chromium 

136 
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Figure 1: Computational Geometry. The prominent GCR incident ions are H, He, C, O, 
Mg, Si and Fe. 
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Figure 2: Boundary Conditions for Aluminum shield. 
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Figure 3: Dose (upper panel) and Dose Equivalent (lower panel) versus Depth in Water 
for Hydrogen on Aluminum shield. 
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Figure 4: Same as Figure 3, except for Helium. 
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Dose Equivalent (cSv/day) Dose (cGy/day) 




Figure 5: Same as Figure 3, except for Carbon. 
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Figure 6: Same as Figure 3, except for Oxygen. 
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Dose Equivalent (cSv/day) Dose (cGy/day) 




Figure 7: Same as Figure 3, except for Magnesium. 
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Figure 8: Same as Figure 3, except for Silicon. 
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Figure 9: Same as Figure 3, except for Iron. 
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Figure 10: Quality factor versus LET. 
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Figure 11: Forward Neutron flux for Hydrogen on Aluminum shield. 
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Figure 12: Backward Neutron flux for Hydrogen on Aluminum shield. 
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Figure 13: Total Neutron flux for Hydrogen on Aluminum shield. 
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Figure 14: Proton flux for Hydrogen on Aluminum shield. 
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Figure 15: 2 H flux for Hydrogen on Aluminum shield. 
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Figure 16: Forward Neutron flux for Helium on Aluminum shield. 
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Figure 17: Backward Neutron flux for Helium on Aluminum shield. 
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Figure 18: Total Neutron flux for Helium on Aluminum shield. 
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Figure 19: Proton flux for Helium on Aluminum shield. 
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Figure 20: 2 H flux for Helium on Aluminum shield. 
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Figure 21: 4 He flux for Helium on Aluminum shield. 
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Figure 22: Forward Neutron flux for Carbon on Aluminum shield. 
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Figure 23: Backward Neutron flux for Carbon on Aluminum shield. 
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Figure 24: Total Neutron flux for Carbon on Aluminum shield. 
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Figure 25: Proton flux for Carbon on Aluminum shield. 
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Figure 26: 2 H flux for Carbon on Aluminum shield. 


42 





H Flux (# particles/cm 2 MeV day) 3 H Flux (# particles/cm 2 MeV day) 


Depth in Water 0 g/crn^ 


Depth in Water 1 0 g/cm^ 



Energy (MeV) 


Energy (MeV) 


Depth in Water 20 g/cm^ 


Depth in Water 30 g/cm^ 



Energy (MeV) 


Energy (MeV) 


Figure 27: 3 H flux for Carbon on Aluminum shield. 
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Figure 28: 3 He flux for Carbon on Aluminum shield. 
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Figure 29: 4 He flux for Carbon on Aluminum shield. 
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Figure 30: Forward Neutron flux for Oxygen on Aluminum shield. 
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Figure 31: Backward Neutron flux for Oxygen on Aluminum shield. 
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Figure 32 : Total Neutron flux for Oxygen on Aluminum shield. 
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Figure 33: Proton flux for Oxygen on Aluminum shield. 
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Figure 34: 2 H flux for Oxygen on Aluminum shield. 
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Figure 35: 3 H flux for Oxygen on Aluminum shield. 
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Figure 36: 3 He flux for Oxygen on Aluminum shield. 
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Figure 37: 4 He flux for Oxygen on Aluminum shield. 
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Figure 38: Forward Neutron flux for Magnesium on Aluminum shield. 
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Figure 39: Backward Neutron flux for Magnesium on Aluminum shield. 
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Figure 40: Total Neutron flux for Magnesium on Aluminum shield. 
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Figure 41: Proton flux for Magnesium on Aluminum shield. 
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Figure 42: 2 H flux for Magnesium on Aluminum shield. 
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Figure 43: 3 H flux for Magnesium on Aluminum shield. 
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Figure 44: 3 He flux for Magnesium on Aluminum shield. 
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Figure 45: 4 He flux for Magnesium on Aluminum shield. 
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Figure 46: Forward Neutron flux for Silicon on Aluminum shield. 
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Figure 47: Backward Neutron flux for Silicon on Aluminum shield. 
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Figure 48: Total Neutron flux for Silicon on Aluminum shield. 
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Figure 49: Proton flux for Silicon on Aluminum shield. 
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Figure 50: 2 H flux for Silicon on Aluminum shield. 
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Figure 51: 3 H flux for Silicon on Aluminum shield. 
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Figure 52 : 3 He flux for Silicon on Aluminum shield. 
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Figure 53: 4 He flux for Silicon on Aluminum shield. 
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Figure 54: Forward Neutron flux for Iron on Aluminum shield. 
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Figure 55: Backward Neutron flux for Iron on Aluminum shield. 
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Figure 56: Total Neutron flux for Iron on Aluminum shield. 


72 




Proton Flux (#particles/cm 2 MeV day) Proton Flux (#particles/cm 2 MeV day) 


Depth in Water 0 g/cm^ 


Depth in Water 10 g/cm^ 




Energy (MeV) 


Energy (MeV) 


Depth in Water 20 g/cm^ 


Depth in Water 30 g/cm^ 




Energy (MeV) 


Energy (MeV) 


Figure 57: Proton flux for Iron on Aluminum shield. 
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Figure 58: 2 H flux for Iron on Aluminum shield. 
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Figure 59: 3 H flux for Iron on Aluminum shield. 
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Figure 60: 3 He flux for Iron on Aluminum shield. 
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Figure 61: 4 He flux for Iron on Aluminum shield. 
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Figure 62: Magnified View of Figure 18 for 30 g/cm 2 depth in water. 
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Figure 63: Magnified View of Figure 20 for 30 g/cm 2 depth in water. 
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Figure 64: Magnified View of Figure 56 for 0 g/cm 2 depth in water. 
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Figure 65: Magnified View of Figure 60 for 30 g/cm 2 depth in water. 
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Figure 66: Heavy ion flux associated with Carbon on Aluminum shield. 
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Figure 67: Heavy ion flux associated with Oxygen on Aluminum shield. 
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Figure 68: Heavy ion flux associated with Magnesium on Aluminum shield. 
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Figure 69: Heavy ion flux associated with Iron on Aluminum shield. 
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and the forward, backward and total neutron fluxes at various depths in the water slab. Comparisons of the secondary light ion 
2H, 3H, 3He and 4He fluxes are also examined. 
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